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Abstract 

In a recent work [l], Anderson et al. used the renormalization group (RG) evolution of the mo- 
mentum distribution to show that, under appropriate conditions, operator expectation values exhibit 
factorization in the two-nucleon system. Factorization is useful because it provides a clean sepa- 
ration of long- and short-distance physics, and suggests a possible interpretation of the universal 
high-momentum dependence and scaling behavior found in nuclear momentum distributions. In 
the present work, we use simple decoupling and scale-separation arguments to extend the results of 
Ref. [1] to arbitrary low-energy A-body states. Using methods that are reminiscent of the operator 
product expansion (OPE) in quantum field theory, we find that the high-momentum tails of momen- 
tum distributions and static structure factors factorize into the product of a universal function of 
momentum that is fixed by two-body physics, and a state-dependent matrix element that is the same 
for both and is sensitive only to low-momentum structure of the many-body state. As a check, we 
apply our factorization relations to two well-studied systems, the unitary Fermi gas and the electron 
gas, and reproduce known expressions for the high-momentum tails of each. 
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I. INTRODUCTION 



Renormalization group (RG) methods play an important role in ab initio nuclear theory 
by extending the range of many computational methods and improving their convergence 
patterns [2HI]. There are numerous RG methods that have been successfully applied to 
nuclear few- and many-body systems in recent years [3J. While the details differ, all such 
methods decouple low- and high-momentum degrees of freedom in a manner that leaves 
low-energy observables invariant. In this paper, we will denote the momentum scale at 
which this decoupling occurs by A. In methods such as the Lee-Suzuki-Okubo similarity 
transformation method or the related V\ ow k approach, A is a floating cutoff beyond which 
high- momentum states have been integrated out [5H7]. For other methods, such as the 
similarity renormalization group (SRG) approach, A gives a measure of how band-diagonal 
the Hamiltonian is in momentum space |8j. In all cases, A serves as a "resolution scale" since 
dynamics above and below this scale are effectively decoupled |3j, E]. 

We emphasize that while observable quantities (such as cross sections) do not change, 
the physics interpretation can (and generally does) change with resolution. It is a com- 
mon misconception that at low-resolution, one is unable to describe phenomena that, at 
high-resolution, are associated with the high-momentum components of low-energy wave 
functions. A prototypical example is the (e, e'p) process at large momentum transfers, where 
theoretical analyses relate such experiments to nuclear momentum distributions if the impulse 
approximation is assumed valid for a high-cutoff interaction [10]. Calculations find nearly 
universal scaling of the high-momentum tails, which is interpreted in terms of short-range 
correlations in the nuclear wave functions Naively it might be thought that this physics 
is beyond the reach of low-momentum approaches, for which wave functions have drastically 
reduced short-range correlations. However, this is not the case: the experimental cross sec- 
tion is unchanged if the corresponding operator is consistently evolved under the RG, even if 
the evolved wave function has almost no high-momentum strength. The formal relationship 
between an operator A ° at an initial high-resolution scale Ao, and the consistently-evolved 
effective operator A at the low-resolution scale A is defined by 



In the (e, e'p) example, one might worry that the consistent evolution embodied by Eq. [T] 
is computationally intractable because the evolved momentum occupation operator might be 
too complicated in practice (e.g., strong non-localities and sizable many-body components, 
etc.). In Ref. [T], some of these questions were addressed by examining the consistent SRG 
evolution of various operators, including momentum distributions and electromagnetic form 
factors in the deuteron. There, all operators were found to flow to smooth, low-momentum 
forms, exhibiting many of the same simplifications as RG-evolved interactions. More inter- 
estingly, under certain kinematic conditions it was found that operator expectation values 
exhibit factorization, which provides a clean separation of long- and short-distance physics 
and an alternative interpretation of the universal high-momentum dependence and scaling 
behavior [T]. 

The proof of factorization presented here and in Ref. [T] follows straightforwardly from 
decoupling and the separation of scales, and is reminiscent of the operator product expan- 
sion (OPE) in quantum field theory. The OPE was developed for the evaluation of singular 
products of local field operators at small separation [T2| IT3] . The utility of the OPE rests 
on factorization; short-distance details decouple from long-distance dynamics. Factorization 



(^°|0 Ao |^) = (^|0 A |^). 
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enables one, for example, to separate the momentum and distance scales in hard-scattering 
processes in terms of perturbative QCD and parton distribution functions. While the meth- 
ods used in the present paper share several similarities with the OPE, a precise connection 
has not yet been made. One key difference is that, in the framework of a local quantum 
field theory, the OPE gives a controlled expansion since the dependence of the Wilson co- 
efficients on the separation r is fixed by the scaling dimensions of the corresponding local 
operators. In the present paper, however, we work in the general domain of non-relativistic 
quantum mechanics (i.e., no assumption of a local QFT). Therefore, we cannot make precise 
statements about the scaling behavior of terms when we expand Fock space operators at one 
resolution scale A in terms of the corresponding operators at another scale A > A. While 
the factorization formulas in the present context are not as controlled as those derived in a 
local quantum field theory using the OPE, they nevertheless provide tools that let us param- 
eterize the high-momentum components of operators which would normally require degrees 
of freedom we do not retain. We can, for example, build effective few-body operators contain- 
ing state-independent functions of high momenta that can be measured directly in few-body 
experiments. These operators can then be employed to make predictions for A-body systems. 

In this paper, we generalize previous developments |T] to derive scaling relations for the 
high-momentum tails of momentum distributions and static structure factors in arbitrary 
low-energy A-body states. In both instances, we find that the expectation value of the 
corresponding operator factorizes into the product of a universal function associated with 
high-momentum (short-distance) physics, and a state-dependent number associated with 
low-momentum (long-distance) structure. The outline of the rest of this paper is as follows: 
In Section [TlJ we review the proof from Ref. jl] that expectation values of high-momentum 
probes factorize in the A = 2 system. In Section |III[ we recast the discussion of factorization 
in a second-quantized language and use it to derive universal scaling relations for momentum 
distributions and static structure factors in general A-body systems. As a test of these 
relations, in Section IV we apply them to two well-studied many-body systems - the unitary 
Fermi gas and the electron gas - to reproduce known expressions for the asymptotic tails of 
the momentum distributions and static structure factors of each system. Our conclusions are 
summarized in Section |Vl and several technical details are relegated to the Appendices. 



II. FACTORIZATION IN THE TWO-BODY SYSTEM 



In Ref. [Ij, Anderson et al. applied renormalization group methods to the two-body 
problem to show that high-momentum components of low-energy wave functions factorize into 
the product of a state-independent function of momentum, and a state-dependent number 
that is sensitive only to low-momentum physics. Using this wave function factorization, it 
is straightforward to show that expectation values of operators that probe high-momentum 
modes similarly factorize into a state-independent piece that encodes the high-momentum 
physics and depends on the particular operator, and a state-dependent number that depends 
only on the low-momentum structure of the state and is identical for all high-momentum 
operators [TJ. As we will show in Section III, the factorization formulas of Ref. [TJ generalize 
to arbitrary low-energy A-body systems, allowing us to derive scaling relations for the high- 
momentum tails of momentum distributions and static structure factors. Since the simple 
factorization in the A = 2 system is the starting point to derive analogous relations for 
general low-energy A-body states, we begin by reviewing the salient points from Ref. [TJ. 
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A. Wave function factorization 



Renormalization group transformations simplify nuclear few- and many-body calculations 
by decoupling low- and high-momentum degrees of freedom leaving low-energy observables 
unchanged [2HI]. In Ref. [I], the analysis was done in the context of similarity renormalization 
group (SRG) transformations, where the resolution scale A provides a measure of how band- 
diagonal the evolved interaction is in momentum spacd^ However, for our present analysis 
we do not have to be very specific about the details of the particular RG implementation. 
All we require is that momentum modes above and below A are effectively decoupled by the 
given transformation. In the center-of-mass frame of the two-body system, this implies that 
the low energy states (\E n \ < A 2 ) are localized in the low-momentum subspacd^ 

pa|v£>«Iv£> fiAlv£>«o, (2) 

where the projection operators Va and Q\ are defined as 

Pa= Tt^IpXpI and Qa = lq>(ql • ( 3 ) 



.0 (27T) 3 '"'^ * A J A (27T, 

Starting from the unevolved Schrodinger equation written in block-matrix form 

VaH^Va VaH^QA (V A ip™\ _ (V A ip™\ ( . 

QaH^Va QaH^Qa) \Qa^J a \QA^J ' 1 ' 

we can solve for the high-momentum projection of any eigenstate as 

Q A |C) = (E a - QaH^Qa^QaH^VaVa |C> 

= (E a - QaHooQa) 1 QaVooVa |Va } , (5) 



oo 
a 



where we have used ("Pa) 2 = "Pa., Hoc = T + Voo, and QaTVa = 0. For low-energy states ip, 
such that \E a \ <C Min[|.E'g#g|] ~ A 2 (where Eqhq are the eigenvalues of QaHooQa), we can 
neglect the E a dependence in Eq. [5] 

/»oo p A -i 

^(q) « - / dq' / dp(q\ |q') V^, p) Cfe) , (6) 

where we've introduced the abbreviation <ig = /f^g-- Assuming for simplicity that ip™ is an 
S-wave state and that the potential V^o(q',p) is slowly varying with respect to p compared 
to V>£°(p) in the region p < A and q' ^> A, we can factorize the low- and high-momentum 
physics by expanding 



/ A #V 00 (q',p)^(p) « ^(q'^OIp^o x f A dp^(p) 
Jo Jo 



1 d 



2 



,A 

X / dpp 2 ^(p) + 
P '=o Jo 



(7) 



1 Using the SRG transformation of Ref. [lj, the evolved potential goes as yj^;^ ~ exp ^— ^ ^ ^ ^ . 

2 For SRG transformations, the high-momentum components of the evolved wave functions are exponentially 
suppressed as exp(— q 4 /A 4 ). The decoupling is exact for RG transformations employing a sharp cutoff A. 
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which gives 

pA pA 

C(q)« 7(q;A) / dp^(p) + rj(q;A) / dpp 2 ^(p) + ... , (8) 
Jo Jo 

where the state-independent functions that carry the q-dependence are defined as 

poo 1 

7 (q; A) = - / dq' (q| |q'> ^(q', 0) , (9) 

1 f°° 1 r/ 2 

i7(q;A) = --/ <ff (q| M^Kc^p') ■ (10) 

2 Ja QaHooQa dp 2 p , =0 

It is known empirically pH [2] that the low-momentum projections of the low-energy 
eigenstates of the bare and evolved Hamiltonians are related by a wave function renormal- 
ization factor Va |Va°) ~ %a |"0a); which reflects the fact that RG evolution does not modify 
long-distance physics. Using that 

/ dpC(p) « Z A / dp^(p) = ^A^W| r=0 = Z A ^(0) (11) 

Jo Jo 
dpp 2 V« (p) « Z A / dpp 2 V£(p) = - Z A V 2 ^(r)| r=0 = -Z A V 2 ^(0) , (12) 



we obtain the momentum space version of Lepage's non-relativistic operator product expan- 
sion [14j relating the short-distance structure of the unevolved or "bare" wave functions to 
those of the low-energy effective theory 

0~(q) « 7 (q; A)Z A ^(0) - 7/(q; A)Z A V 2 ^(0) + • • • . (13) 
If we keep only the leading term in the expansion 

C(q) « 7(q; a)z a ^(o) = 7 ( q; a)z a tdp^ip) , (14) 

Jo 

we see that the high-momentum components of the low-energy eigenstates are factorized into 
a state-independent function 7(q;A), which summarizes the short-distance behavior of the 
wave function, and a state-dependent coefficient that probes the low-momentum structure of 
the state. 



B. Effective operators and factorization 

Given the wave function factorization in Eq. [I4j we can now derive analogous factorization 
formulas for expectation values of general operators in the A = 2 system. Consider the 
expectation value of an operator O in a low-energy eigenstate of the unevolved Hamiltonian 

pA pA pA pod 

W\6K)= / dp ^ / vr(p)o(p,P / )^(p') + / dp «r(p)0(p,q)dq) 

Jo Jo Jo Ja 

poo pA poo poo 

+ dq dpV>r(q)0(q,p)«p) + / ^ dq' VE* (q)O(q, q')^~(q') , 
Ja Jo Ja Ja 



(15) 
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where we have explicitly separated the low- and high-momentum integrals in forming the 
matrix element. Next, we insert Eq. 14 and ^£°(p) ~ ^A.V'a(p) f° r ^ ne high- and low- 



momentum components of ip^, respectively. Since the matrix elements 0(p, q) and 0(q, p) 
involve well-separated momenta, we perform a Taylor expansion about p = and keep only 
the leading term giving 



(CMC)«^i [ A dP /V^*(p)c>(p,p')^(pO 

Jo Jo 

POD 

+ 2Z 2 A \^(0)\ 2 / dgfO(0,q) 7 (q;A) 



A 

oo 



POD POD 

+ Zl\^(0)\ 2 dq rfg' 7 *(q;A)0(q,q')7(q , ;A)- 
J A J A 

(16) 

Since the evolved wave functions ^(k) have vanishing or exponentially suppressed support 
for k > A, we can re- write this as 

md\1>?) « Zl(^\d\^) + <7 (0) (A) (^l^(r)l^) , (17) 
where the coupling g(°\A) is defined as 

POO 

0<°>(A)=22fl / dgO(0,q) 7 (q;A) 
JA 

+ Zi / dg/ rfg' 7 *(q;A)0(q,q') 7 (q / ;A). (18) 
Ja Ja 

Recalling that the consistently evolved effective operator is defined by 

(CIOIC) = «|o A |^), (19) 



we see from Eq. [17] that 

A « ZiO + ^(A)^(r) + ... , (20) 
where the ". . ." contains higher derivatives of delta functions that arise from the gradient 



terms in Eq. 13 as well as higher-order terms in the expansion of 0(q, p) about p = 0. 
In this way, we see that the RG-evolved operators take on a universal form; the effects 
of the integrated-out high-momentum modes are absorbed in a rescaling of the unevolved 
operator at the initial resolution scale, plus a series of local, state-independent corrections 
that take the form of a derivative expansion with A-dependent couplings [14, 15J. As stressed 
by Lepage |14| . the universal form of these local corrections is analogous to the multipole 
expansion in classical electromagnetism; just as multipole moments may be calculated from 
an underlying theory (e.g., the true charge and current densities) or extracted from a finite 
number of experimental data, the same holds true for the couplings Z\, g(°'(A), etc. 



Let us now consider the implications of Eq. 17 for operators that predominantly probe 



high-momentum components of low-energy states. Since such operators have negligible 



strength at low-momentum V\OV\ ~ 0, the first term in Eq. 17 vanishes, leaving 



(ClOlC) ~ S (0) (A) (^ {3) (r)l^> • (21) 



6 



Therefore, the expectation value of any operator that probes the high-momentum struc- 
ture of low-energy states factorizes into a state-independent piece, g(°\A), that depends 



on the particular high-momentum operator via Eq. 18, times a state-dependent number, 
("0a I ^ WlV'a)) that is the same for any high-momentum O, and is only sensitive to the 
low-momentum structure of the state since T^aIVq ) ~ Z^\ip A ). 

The momentum distribution n q = a q a q for q ^> A is a prototypical example of an operator 
that is sensitive to the high-momentum structure of wave functions. Since n q = |q)(q| for 



the A = 2 system in the center-of-mass frame, Eq. 21 becomes 



« r (q; A) Zi m0)\\ (22) 

We see that momentum distributions in all low-energy states (\E a \ < A 2 ) in the A = 2 system 
share the same q-dependence for q > A. In fact, we will find in the following section that 



the factorization formula, Eq. 22 , generalizes to arbitrary A-body systems. 



III. FACTORIZATION IN THE ,4-BODY SYSTEM 
A. Evolved creation and annihilation operators 

In order to proceed beyond the A = 2 system, it is convenient to recast the results from the 
previous section in a second-quantized language. As a first step, we examine how the Fock 
space creation and annihilation operators evolve under RG transformations. Suppressing 
non-essential spin and isospin indices, the transformed operators can be expanded on the 
original operator basis as 



a q ^ — a q 



+ O k i> k 2 )aj Ci a| C2 a kl+k2 _ q + . . . = a q + 5a^ (23) 

ki,k 2 



where the ". . ." contains higher- rank terms that are generated (a^a^a^aa, etc.) when the RG 
evolution is carried out beyond the two-body level, i.e., when induced 3- and higher-body 
interactions in H A are not truncated during the flow. Note that the form of the coupling 



function C q v (ki,k2) can be constrained further by boosting both sides of Eq. 23 and using 
Galilean invariance to write 



q-P - a q-P + Cq( k l' k 2)a kl _ P a k2 _ P ak 1 +k 2 -q-P 

ki,k2 

= a q-P + C i ^ + P > k 2 + P)a kl a k 2 a k 1+ k 2 -q+P 

ki,k 2 

= 4-P + C 'q V -P( k l' k 2) a k 1 a k 2 a k 1 +k 2 -q+P , (24) 



ki,k 2 

which implies 

C^ki+P^ + P) = ^(ki.ka). (25) 

In the following, we restrict our attention to the leading non-trivial term in Eq. [23j This 
corresponds to neglecting induced three- and higher-body interactions in H A since the co- 
efficient function Cq(ki,k 2 ) is uniquely determined from the RG evolution in the two-body 
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system. This can be seen by considering the following matrix element between the zero- 
particle vacuum (which does not evolve under the RG) and a two-body eigenstate in the bare 
and evolved theories 

(Cl4 +P 4 10) = (^iai^a^JO) 

2 ~r 2 y 2 ' p 2 ' p 



<v> T p +P 4-pi°> + «i^ T i + p aT f -pi°> 



A la f 

a\ u I r 
Ait 



f +p a f-p 



P, ^ P 



^4, 4 
2 ^ p 2 

|0> + E ^^(P/ 2 + k, P/2 - k) |4 

k 



, at , 

+k ^-k 



10) 
(26) 



where <5a^|0) = was used in the second line and Eq. 25 was used in the last step. Since the 
dependence on the COM momentum P cancels on both sides, we are left with 



^r(p) = v£*(p) + E - k ) ^*( k ) . 

k 

which can be inverted using the completeness of the {ip^} t° giv^] 

c p A (k,-k) = ^(k|i A )(cip)-5 k , P . 



(27) 



(28) 



One can use the results of Section[TT]to evaluate two important limiting cases of Eq. [28] that 
will prove useful below. First, consider C A (p', — p') for p,p' < A. Using VaI^™) ~ Z\\x/j a ) 
for low-energy states and PaI^q) ~ for \E a \ > A 2 , Eq. 28 becomes 



Cp A ( P \-p') E (p'IV&^Ip)-*, 

l^l<A 2 

« (Z A -1) V,p- 



PiP 



(29) 



In the last step, we used that the low-energy evolved eigenstates span the low-momentum 
subspace due to decoupling 

Pa = X>><pI w E IV£XV£|- (so) 

P<A l^ Q |<A2 

The other important limiting case is C A (p, — p) for p < A and g > A. Inserting Eq. 
Eq. 28 then gives 

C A (p,-p)«Z A7 (g;A) E (P'l^><^|r = 0) 

\E a \<A- 2 

nZ Al (q;A), (31) 
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into 



where Eq. 30 was used in the final step. 



3 It is helpful to think of the SRG or the Lee-Suzuki-Okubo similarity transformation method, where the 
decoupling of low- and high-momentum modes is accomplished by a unitary transformation. By unitarity, 
one has Yl IV^XV^I =1 = 2 I'^aXV'a l> where the sum over a is unrestricted. 

a a 
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B. Factorization for momentum distributions 



In Eq. |22l we found that the expectation value of the momentum distribution in low- 



energy two-body states factorizes for large q > A. Using the second-quantized formulation 



of Section HI A we will now show that a similar factorization occurs for general low-energy 
A-bodj states. We begin by considering the expectation value of the consistently evolved 
momentum distribution operator for q ^> A in an arbitrary low-energy A-body state 



= + K A)t « q + 4M A) + H A) H A) }i<J • 



(32) 



This is an exact equality provided that a) the evolved Hamiltonian H A includes all induced 
3-, 4-, . . . v4-body interactions generated by the RG evolution, b) all higher-order terms for 
<W A)t and <W A) in Eq. " 



23 



are included, and c) all possible 1-, 2-, . . ., A-body operators 
generated by the terms in the curly brackets are keptQ However, since we are only interested 
in the high-momentum tail of Eq. 32 and since one expects induced 3-body and higher 
operators contributing to H A and [OqtJq]^ to be subleading so long as one doesn't evolve 
too low in A [3], we will neglect them. 

In what follows, we assume A is of the same order as the physical momentum scales 
that characterize A ( e -§-' ^ e F errm momentum, hp, for homogenous systems, ^Jmuj/h 
for harmonically-trapped systems, etc.). Due to decoupling, the low-energy states ipa A have 



vanishingly small support at high momentum. Therefore, any term in [ala q ]^ that annihi- 



lates a high-momentum particle from \ipt A ) or (^A 
find that Eq. 32 becomes 



A 

Ot,A I 



will be suppressed. In this limit, we 



n q « 



<<jH A)t K A) i<A 

K 

2" 



A 



A 



2. c « 

k,k',K,K' 



+ k, 



) 
K 

~2 



k\Ci 



K' 

Y 



k' 



K' 



x (V> 



a, A i a K +k a K_ k a K-q a K'-q a ^- 



k' 



£ c i (I + k | - k ) c i (I + k/ > I - k ') «J4 +k 4_ k « f ^f-,iO 

:.k'.K ^ - / V J - - 

£ ( k ' - k ) C t* ^ " k ') (<j4 +k 4_ k «f +k^f-k'l< A ) , (33) 



k,k'.K 



where we have anti-commuted ctK-q to the right and dropped the normal-ordered three-body 
term in going from the second to third line. The low-momentum nature of ip A implies that 
dominant terms in the sum are for |K/2 ± k| and |K/2 ± k'| < A. Consequently, we have a 
mismatch of scales |q — K/2| S> |k|, |k'|, which together with Eq. 31 gives 



7 2 (q-K/2;A)(^Ja t f+k a t f _ k aK +k ,aK_ k ,|< A ) 

k.k'.K 

~ Zfo 2 (q; A) (<Al a k +k 4-k a f+k' a f-k'l<J , 



(34) 



k,k',K 



4 An A-body operator is defined as a normal-ordered string of A a^'s and A a's. 
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where we've used q 3> K/2 in the last step^J In this way, we see the large-g tails of momentum 
distributions for arbitrary low-energy /L-body states share the same universal g-dependence. 



In nuclear physics, Eq. 34 provides an alternative to the usual explanations based on short- 



range correlations [TOl [T6] as to why calculated momentum distributions in various nuclei 



and nuclear matter scale with each other at large q. In Section IV, we will use Eq. 34 and the 



analogous expression, Eq. |4TJ to reproduce known asymptotic expressions for the momentum 
distributions and static structure factors for two well-studied many-body systems, the unitary 
Fermi gas and the electron gas. 

C. Factorization for static structure factors 

The static structure factor is an important quantity that contains information about 
density-density correlations in a many-body system. For a many-body system of fermions 
with two spin states, the correlations between the densities of the two spin states are par- 
ticularly important. The corresponding static structure factor S^(q) for a homogeneous 
system is the Fourier transform in the relative coordinate r x — r 2 of the density correlator 
(^a' A \P-\{ Y ^)Pi{ Y 'i)\' l Pa' A )- Using similar arguments as for the momentum distribution, we now 
show that at large momentum 5^(q) factorizes into a universal function of q times a matrix 
element of a delta function in the evolved low-momentum wave functions. 

Starting from the definition of S^(q) in the unevolved theory 

■%(q) = (c;>l(q)a(q)ic; A > = E^J a p,t a p+ q ^ +q ,,v,J^) 

p,p' 



P,P 



we consider the expectation value of the consistently evolved operator S^(q; A) in the evolved 



wave functions. Using Eq. 23 for the evolved creation/annihilation operators, we have 



(CJ%(q)ICJ = (<J%(q;A)|< A > = (<J(%(q) + ^(q))l<,>- (36) 

This is an exact relation only if all induced many-body operators (up to rank-^4 for the A-body 
system) are kept in H A and 5S A (q). As with our analysis of the momentum distribution, we 
neglect these many-body contributions by a) restricting the expansion of a*- A ^ and a^ A ' to 



the leading terms shown in Eq. 23 and, b) truncating 5S A (q) to two-body operators 
<^(q)~ Yl C q A +k'( k > " k ) a K +k /K_ k /K +k , 4 aK_ k , t + h.c. 

K.k.k' 

+ Yl ^P+q-f ( k >- k )C£-f ( k '> - k ') a |+k/f -k/f +k',| a f "k'.t 
P.K.k.k' 

= 8S A (q) + 5S A (q) , (37) 
where 5S A ^(q) denotes the terms linear (quadratic) in the expansion coefficients C A . 



5 For systems where 7(q; A) exhibits a power-law decay, the corrections for non-zero K do not modify the 
power-law tail of n q . 
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With the approximate form of the evolved operator in hand, we can now evaluate Eq. 36 



for q A, where once again A is assumed to be of the same order as the physical scales that 
characterize the system. The expectation value of the bare operator, S^(q), in the evolved 
low-momentum wave functions is negligible since it involves the removal of a high-momentum 
particle. Therefore, we have 

(C;j%(q)IC> « (<J(^(q) + ^2(q))l< A ) 

= 2 E ^+k'( k '" k ) (^ A |a t K +kit «k_ k ^K +k/ji aK_ k;it |^ A ) 



K,k,k' 2 



P,K,k,k' 

(38) 

Due to the low-momentum structure of the evolved wave functions, it is clear that the sums 
over momenta K, k, k' are effectively cutoff at A, while the summation over P is unrestricted 
in the second term of Eq. 38 Performing a Taylor series expansion of the coefficient functions 
in powers of the small momenta K, k, k' and keeping just the leading term gives 



<V£j3&(q)llC> w 1 2 ^ A ( ' ) + E C P +q (0. 0)^(0,0) 

E (^j^f+k/f-k/f+k'^f-k'.ti^J 



p 

X 



K.k,k' 



(39) 

To proceed further, we consider the following three regions that arise in the sum over P: 

• Region I): |P + q| > A and |P| > A 

• Region II): |P + q| < A and |P| > A 

• Region III): |P + q| > A and |P| < A. 

Regions II) and III) are trivial since the C A coefficients involving all soft momenta give a 
delta function, Eq. [29} that allows the sums to be performed. Together with Eq. |3lJ we have 

^C£ +q (0,0)C£(0,0) ^^(Z A -l)5p, q Z A 7(P;A) = Z A (Z A - l) 7 (q; A) 
p,ii p,ii 

^c^ +q (o,o)^(o,o) ^^(z A -i)5p, ZA7(P + q;A) = z A (z A - i) 7 ( q; A) , 

P,III P,III 

(40) 

which gives 

(CJ^(q) IO « I 2 ^7(q; a) + E z l + * A )^( p - A ) 



PJ 

X 



E (<J a K +k , t a f -k,/f + k',^f-k', t i<.) • 

K,k,k' 

(41) 

As with the momentum distribution, we see that the high-momentum tail of the static struc- 
ture factor in a general low-energy v4-body state factorizes into a universal function of q, 
multiplied by a state-dependent matrix element that is controlled entirely by low-momentum 
physics. 
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D. Factorization for general high-momentum operators 



While our explicit proofs of factorization have thus far been limited to the momentum 
distribution and the static structure factor, the phenomena is very general and can be qual- 
itatively understood from Eq. [T] Consider an operator at the initial high-resolution scale Ao 
that probes high-momentum modes, Oq°, where the subscript q indicates that the second- 
quantized expression involves creation (annihilation) operators that add (remove) a high- 
momentum particle. We assume q is much larger than any physical scale that characterizes 
the low-energy state ipn°i an d we also assume |q| <C A so that the expectation value of 
Oq° is non-vanishing. Now consider the consistently evolved Oq, where A |q| Ao, and 
expand it as a polynomial in creation/annihilation operators defined at Ao- Schematically, 
we have 

6^ = J2^A a (42) 

a 

where A a denotes a normal-ordered string of creation/annihilation operators at Ao, a is a 
collective index for the different momentum modes being added/removed, and g® is a c- 
number coefficient. Inserting this into Eq. [T], we have 

(^°l^°|^ ) = ^<(^|i a |^). (43) 

a 

Due to the low-momentum nature of the evolved wave functions, we find that only the A a 



involving the addition/removal of low-momentum (< A) modes contributes to Eq. 43 Since 
all momentum modes contained in a obey k a /q <C 1, the c-number coefficients can be Taylor- 
expanded in the soft-momenta. Note that this expansion should be well-defined since the 
#q encode contributions from loop integrals that have both ultraviolet (Ao) and infrared (A) 
cutoffs in place, thus preventing any singular behavior from arising. In this way, we see that 
the universal q-dependence factorizes, and the remaining state-dependence is given by matrix 
elements of low-momentum operators. 



IV. EXAMPLES 



As a check of our factorization formulas Eq. [34] and Eq. [4T| we use them to reproduce 
known expressions for the high-momentum tails of n q and ^(q) for two well-studied systems, 
the unitary Fermi gas (UFG) and the electron gas. 



A. Unitary Fermi gas 

1. Momentum distribution 

In the case of a unitary Fermi gas described by a contact interaction, the coefficient 7(q; A), 
and hence the large momentum tails of n q and S(q), can be calculated analytically. Consider 
the two-body Hamiltonian with a spin-independent contact interaction 
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where A is the ultraviolet cutoff on all momenta of the theory. Here, we assume that A is 
much larger than any relevant low-energy scales in the problem such as the inverse scattering 
length or the Fermi momentum. The coupling constant g(A Q ) is determined by matching the 
scattering amplitude at threshold to the S'-wave scattering length a and is given by p2] 



0(A„) 



An 



Ana 2tt 2 



-i 



(45) 



To obtain an explicit expression for 7(q;A) in Eq. [9j the operator (Qa-^ooQa) 1 can be 
constructed with the aid of the operator identity 



1 



1 - A- l B + A- X BA- X B - . . OA" 1 



A + B 

where A — > QaTQ a and B — > QaVsQa giving 

(q|(Q A ^ooQA) |q) 



(46) 



7(q;A) 



2m 

g(Ap) 
q 2 



(27T) 



! _£(Ao) [ A0 <ti4_ + 



q" 1 J a 27r 2 g' 2 
2tt 2 



<?(Ao) 



E 

ra=0 



g(Ao) ■ (Ap - A) 
2vr 2 



g 2 27T 2 + ^(A )-(A -A)- 



This can be simplified further since Eq. [45] implies 



27r 2 + s(A )-(A -A) 
which gives 



Ap-A 
2tt 2 



7(q;A) 



1 


-i 


" 1 A 


-l 


<?(A ). 




47ra 27r 2 





<?(A) 



ff(A) 



Inserting into Eq. |34j we find 

9 



^A5 ,2 (A) / /A |t t I / A \ 

— a (^, A l«K +k «^_ k a f+k/ aK_ k ,|^) 



(47) 
(48) 

(49) 
(50) 



k,k'.K 



where A is of the same order of magnitude as the relevant low-energy scales of the system 
and A « g « A . 

In Ref. |17| . Braaten and Platter used the operator product expansion to show that the 
tail of the UFG momentum distribution behaves lik^3 



Tlr 



C(A C 



(51) 



k,k',K 



where C(Ao) is often known as Tan's contact parameter. In the Appendix, we will show that, 
at the level of approximating the evolved creation/annihilation operators by the leading-order 



6 Shina Tan provided the first derivation of Eq. 5f using generalized functions |18| . Many different derivations 



can be found in the literature, see Ref. 1 19 for details. 
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expression in Eq. 23 and truncating induced three- and higher-body operators, the following 
relationship holds 



Z 2 A C(A) = C(Ao) 



(52) 



Therefore, Eqs. 50 and 51 are equivalent at the level of approximations made thus far. 
Heuristically, we can understand this equivalence since we expect that — > 1 as A — > Aq. 



2. Static structure factor 



Turning next to the asymptotic expression for the static structure factor, Eq. 41 our task 
is to evaluate the following term 



5>i T (P + q;A) 7 (P;A) Z\g\A) f 



d 3 P 0(|P 




-A)0(|P + q 


-A) 


(2tt) 3 


P + qp 


! |p 


2 



(53) 



where we've taken the infinite volume limit to convert the sum to an integral and substituted 



Eq. [49] for 7. Note that the integral in Eq. [53] has an implicit ultra-violet cutoff Ao ^> A. 
To evaluate the integral, we use that there are two regions for which the theta function 
#(|P + q| — A) = 1 independent of P • q, and one region where it depends on angle to write 



d 3 P 6{\P 




-A)0(|P 


+ 


q 


-A) 


(2tt) 3 


P + q 


2 


P 


2 



'high 



medium 



+ I] 



low 



(54) 



For |P| > |q|, we have |P + q| > |P| — |q|, which implies that for |P| > A + |q|, then 
0(\P + q| — A) = 1 independent of P ■ q. In this case, the limits of the angular integration 
are unrestricted 



'high 



1 

4^ 



A rl 

dP / dx 

A+q J-l 



1 



Aq— ^OO 



1 



P 2 + q 2 + 2Pqx 



An 2 q 



Lig 



Q 



q + X 



- Li 2 



q + X 



(55) 



where Li2(x) is the polylogarithm function. Using that A/q ^ 1 and keeping just the leading 
term gives 



high 



1 


V 


A, 

h - log 




A 


A-K 2 q 


_T " 


(Q 


< l. 



(56) 



Similarly, for |P| < |q| the limits of the angular integration are unrestricted if |P| < |q| — A 
giving 



'low 



47T 2 



47r 2 g 
1 

An 2 q 



1 



q-A rl 

dP / dx — — 
A 7-i P 2 + q 2 + 2Pqx 



2 



Li 2 



A 



7T A, [A 
4 + i bg U 



3A 
Q 



(57) 
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where we've used A/q <C 1 in the last step. Finally, we consider the intermediate region 
q — A < P < q + A where #(|P + q|— A) = 1 places restrictions on the the limits of the 
angular integration. In this case, the theta function requires x > x m i n = 



2Pq 



' medium 



1 
1 

4ir 2 q 
1 

47r 2 g 



dP 



log 

2 A 
. Q 



log 



dx 



1 



log 

2q 
A 



Inserting Eqs. 56 58 in Eq. |54| gives 



P 2 + q 2 + 2Pqx 
l + A/q 



1- A/9 



Li 2 



Li 5 



(5c' 



4n 2 q 



7T 

y 



4A 
Q 



(59) 



which together with Eq. [41] and Eq. 49 yields 

2 



%(q) 



g^(A) 
2 

~g 2 ^(A) 
l l 

8g 



+ 



+ 



1 

8q 
1 

8q 



2-Kaq^ 



A 

TT 2 q 2 

A 
7T 2 q 2 



C(A C 



Z\C{A) 



C(A ) 



(60) 



where we used Eq. 52 and the explicit form of the coupling g(A), Eq. 45, in the second and 
third lines, respectively. As with the momentum distribution, Eq. 60 agrees with the known 
result that has been previously derived by a number of different methods |19j . 



B. Electron gas 

1. Momentum distribution 



As our second check of Eq. 34 and Eq. 41 , we derive the large-momentum limit of the 



momentum distribution and static structure factor for Coulombic systems. Unlike the unitary 
Fermi gas, we were unable to evaluate 7(q; A) in closed form. Therefore, we turn to a 
perturbative calculation and expand the Q-space propagator 



1 -v^ + ^v^v- 1 



QaHQ a Q a TQ a Q a TQ a Q a TQ a Q a TQ a Q a TQ a Q a TQ / 



which together with Eq. [9] gives the first- and second-order contributions to 7 



(61) 



A (2tt) 3 q 2 ^ ^ J q' 2 
Air 



a q 4 



(62) 
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and 



7 (2) (q;A) = - 



d?q' 1 47re 2 47re 2 
(2ir) 3 q 2 q' 2 |q — q'| 2 g /2 



a 2 g 4 A 



(63) 



where we've kept the leading term in 1/q for the second-order contribution and = is 
the Bohr radius. We assume that perturbation theory is justified provided 



7 



(2) 2 1 



7W 7T OqA 



< 1 



A> 



2 1 

7T Oq 



and restrict our attention to the leading term, Eq. 62 Inserting this into Eq. [34] gives 

16tt 2 

k,k'.K 



zi E <<>Ut 



^K-k^f+k' 



,t a f-k'J< A > 



(64) 



(65) 



Apart from the factors and the evolved wave functions ^ , this is very similar to the 
known result first derived by Kimball [20J 



167T 



E <C>f + k,t a W a f+ k '*t a f -k',J^> 



(66) 



k,k',K 



As with the unitary Fermi gas, one can make a heuristic argument that Eq. [65] and Eq. [66 

o 



are equivalent since Za — > 1 and ^ — )■ as A — >■ 00. More precisely, we will show in the 
Appendix that 



k,k',K 



^{i+°(^)}E<<J<W*- 

1 \ 0/ ^ k,k'.K 



k,^f+k 



(67) 



so that Eqs. 



65 



and 



66 



are equivalent up to terms of order 0(-^), which are presumed to 



be small by virtue of~Eq. 64 



2. Static structure factor 



Turning next to the application of Eq. [41] to Coulomb systems, our task is to evaluate the 
following term 

»(P + q ;AMP;A) . (68) 

As before, we split the integral into a sum of three terms 



f d*P fl(|P|-A)fl(|P + q|-A) _ 

J (27r) 3 |P + q| 4 |P| 4 -'high ~r J medium ~r J low V 0J J 
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where T^igh corresponds to |P| > A+ |q|, I me dium corresponds to |q| — A < |P| < |q| + A, and 
I\ ow corresponds |P| < |q| — A. For Thigh an d -^low; the angular integrals are trivial 



'high 



1 

471-2 



P 2 dP 



A+q 



P 4 (P 2 + q 2 + 2PqxY 



4tt 2 2Aq 4 



O 



' low 



1 

1 



qr-A 



4tt 2 2Ag 4 



P 2 dP 



+ O 



-i P A (P 2 + q 2 + 2PqxY 



For L 



medium' 



the Heaviside theta functions restrict the angular integral 



'medium 



1 

Air 2 
1 1 



q+A 



P 2 dP 



q-A 



1 



1 



P A (P 2 + q 2 + 2PqxY 



where x r 



A 2_p2_ q 2 



2Pq 



%(q) 



8%Zl 
a q 4 

8vr 1 
a Q q 4 



4vr 2 Ag 4 
. Combining Eqs. 
2 1 



+ O 



68- 



72 



1 - 



7r Aa c 



with Eq. 



41 



we find 



K,k,k' 



K,k,k' 



(70) 



(71) 



(72) 



(73) 



where we've used Eq. 64 and Eq. 67 to obtain the second line. Once again, this is in agreement 
with the previously known result of Kimball 



V. SUMMARY 

In this paper, we have used elementary RG arguments to show that, for general low- 
energy many-body states, the high-momentum tails of momentum distributions and static 
structure factors factorize into the product of a universal function of momentum that is 
fixed (in leading-order) by two-body physics, and a state-dependent matrix element that 
is sensitive only to low-momentum structure of the many-body state, and is the same for 
both. This generalizes the results of Anderson et ai, who derived analogous relations in the 
two-body system [lj, and suggests a possible interpretation of the universal high-momentum 
dependence and scaling behavior found in nuclear momentum distributions in the analysis of 
(e, e'p) reactions. As a check, we have successfully applied our factorization relations to two 
well-studied systems, the unitary Fermi gas and the electron gas, reproducing known results 
for the high-momentum tails of each. 

Our proof of factorization follows from decoupling and the separation of scales, and re- 
sembles aspects of the OPE in quantum field theory. Unfortunately, we have not been able 
to establish a precise connection. The main difference appears to be that, in a local quantum 
field theoretical framework, the OPE offers a controlled expansion since the scaling dimension 
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of a given local operator uniquely fixes the r-dependence of the corresponding Wilson coeffi- 
cient, making the truncation of the expansion controllable. In contrast, in the present paper 
we work in the domain of general non-relativistic quantum mechanics and do not require that 
the system is described by a local QFT. This relaxation of assumptions allows us to extend 
the notion of factorization and OPE-like methods to a wider class of problems, albeit in a 
less controlled fashion since we cannot make precise statements about the scaling properties 
of the operators kept/omitted in our expansions. Nevertheless, the methods presented in this 
paper may still be useful in low-energy nuclear physics, as they provide tools that allow us 
to parameterize the high-momentum components of operators that would normally require 
degrees of freedom that we do not retain. We can, for example, build effective operators con- 
taining state-independent functions of high momenta that can in principle be extracted from 
few-body data, and subsequently used to make predictions for high-momentum processes in 
A-body systems. 
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Appendix A: Tan's contact at different resolutions 



In this appendix, we derive the following relation (Eq. 52 in the text) 

Z 2 A C(A) = C{Ao) , 
where Tan's contact parameter is defined as 



A 



C(A ) = <7 2 (A ) E (^l^+k^f-^f+k^f-k'l^) 

K,k,k' 



and g(Ao) is given by Eq. 45 Let us begin with the following relation 

Ao Ao 



E ^t\ a l + A-* a f + v a f-^t) = E « 



K,k,k' 



K,k,k' 



a f+k a f-k 



(Al) 



(A2) 



l A 



~2 +k ~2~ 



(A3) 



where the evolved pair creation operators are given (within the same approximations made 
in the text) by 



Ok 



ate 



= a 



E 

p 



Ck (P> -P) 4 , n a k 

2 ~ V 2 " 



r 



+ 6a 



= +k 



Ok 



(A4) 
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and similarly for the pair annihilation operators. Note that the form of Eq. A~4 is dictated by 
the two approximations made in the text, namely i) the neglect of higher-order corrections 
in Eq. 23 for the evolved creation/annihilation operators and ii) the neglect of induced three- 
and higher-body operators in the evolution. At this level of approximation, the right-hand 



side of Eq. |A3| becomes 







Ao r 


RHS of Eq. 


A3 


= E 






K,k,k' . 



(^L I fl f +k fl f -k fl f +k' fl ^ -k' I V>a, A > 



k k 

2 1 2 



ok u k ,aK_ k ,|^ iA ) 

+ « A I «K +k ak K + J&I K _ k , | ^ A ) 



(A5) 



r 



+ (^, A l 5at K +k at K_ k a f+k'^f-k'l^, A ) 
Next, we split each summation into low- and high-momentum regions and use decoupling 



together with the asymptotic forms in Eqs. 29 and 31 to simplify the resulting expressions. As 



a consequence of decoupling, the first term in Eq. A5 stays the same but with the summations 
effectively cutoff at A <C Aq. The second term is given by 



A A 



E 52 C £(P>-P) ^aJ a f +p a f- p a f+k' a f-^aJ 
K.p k' k 



(A6) 



where the sums over K,p, and k' are cutoff at A due to decoupling. Splitting the sum over 
k into low- and high-momentum regions and using Eqs. 29 and 31 for C k (p, — p) gives 



Ao 



Yl (VClMi +k 4<_ k a f+k' a f-k'l^aJ 
K.k.k' 



(A7) 



(Z A - 1 + Z A 7(A)) «J4 

K,k,k' 



a k_ k aK +k ,aK_ k ,|^ Q ) 

2 1 2 ^ z 



where we've defined 



Ao 

7(A) = X)7(q;A). 

q=A 



(A8) 



One easily finds that the third term in Eq. |A5| gives the same contribution as the second 



The fourth and final term in Eq. A5 is given by 

A A 

E E a k(p, -p)Ck A '(p', -po «Ji P i P «i + /f^it) • (aq) 

K,p,p' k,k' 
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As before, decoupling implies the sums over K, p, and p' are cutoff at A. Splitting the 
unrestricted sums over k and k' into low- and high-momentum regions then gives 

E (i A J^_ k ^f-.lt) « ((Za - l) 2 + 2Z A (Z A - 1) 7 (A) + ^7 2 (A)) 



E (^,Al a k +k a k_ k aK +k ,aK_ k ,|^ A ) . 



K,k,k' 



(A10) 



Collecting terms and simplifying, Eq. |A3| becomes 



A„ 



E (<Al4 +k 4-k a f+k' a f-k'l< Q A )~ ^i(7(A) + i) s 



K,k,k' 



(An; 



X 



E A I a l< +k a l< - k fl f +k' a f -k' I ^a, A ) 



K,k,k' 



For the unitary Fermi gas, we can evaluate 7 analytically using Eq. [49 

* A ° d 3 q -g(A) g(A) 



A 



7(A) = ^ 7 (q;A) 

q=A 



— >■ 



a (2vr) 3 g 2 



2tt 2 



(A - A) . 



(A12) 



Inserting this into Eq. All and using the identity 

<?(A) , <?(A) 



2(A C 



2tt 2 



(A - A) 



(A13) 



one finds the desired relation 



A 



C(A )=g 2 (A ) E «Al4 +k a k_ k a f+k'«f-k#a 



A 



K,k,k' 
A 



z Wi A ) E (^Al a k + k a k_ k a f+k' a f-k'l^ A ) 



(A14) 



K.k,k' 



ZlC(A) . 



Appendix B: Coulomb pair-distribution function at different resolutions 



In this appendix, we derive the following relation (Eq. 67 in the text) for the Coulomb 



gas, 



E (^Al G f +k at K-k a f+k' a f-k'ICA) « 

k,k',K 

+ (aL) } S «>Wk4-k a ¥ ^..10 • 

V U/ } k,k',K 



(Bl) 
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Since the proof closely mirrors what was done for the unitary Fermi gas, we begin from 
Eq. |All| with A oo 



E «j4 +k 4- k a f +k'«f -k'ld ~ Z l (7(A) + 1) 



K,k,k' 



(B2) 



K,k,k' 



Using the leading perturbative expression for 7(q; A) in Eq. 62, we have 



7(A) = ^ 7 (q;A) 

q=A 



d 3 q -4tt 1 2 1 



A (27r) 3 a g 4 7r Aa 



(B3) 



Inserting this into Eq. B2 then gives 



°° / 2 1 \ 

E (CJ4 + k4-k a f +k ' a f- k 'io - ^ ( 1 - -xj 

K,k,k' ~ \ u/ 



X 



7T Aao, 

E (^Ai a k + k a i-k a f+k' a f-k'K A ) 



(B4) 



K,k,k' 



which is equivalent to Eq. Bl since Aao ^ 1 by assumption. 
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